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Abstract 



We study, within the Hartree approximation, the effects of the electron- 
electron interaction on the energy spectrum of a two-dimensional electron 
gas in a perpendicular homogeneous magnetic field and a lateral superlat- 
tice potential with square symmetry. Due to the strong screening effects, 
the bandwidth of the Landau bands oscillates strongly with their filling. For 
short enough periods and strong enough modulation of the superlattice po- 
tential, the miniband structure of the Landau bands can be resolved in the 
thermodynamic density of states. 
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I. INTRODUCTION 



The study of the characteristics of a two-dimensional electron gas (2DEG) 
in a homogeneous perpendicular quantizing magnetic field B and a lateral pe- 
riodic superlattice potential has been revived, since due to recent technological 
developments it now is possible to manifacture such systems with controlable 
properties in semiconductor heterostructures. Measurements of the magne- 
totransport in superlattices weakly modulated in one or two directions have 
revealed novel oscillations in the conductivity, the Weiss oscillations, reflect- 
ing the commensurability of the cyclotron radius R c = l 2 kp = vp/uj c of the 
electrons at the Fermi energy Ep = h 2 k F /(2m*) and the modulation period 
L, I = (ch/eB) 1 / 2 is the magnetic length. [j^,|2| The Weiss oscillations are 
superimposed on top of the Shubnikov-de Haas oscillations which are well 
known in homogeneous 2DEG's and caused by the commensurability of the 
Fermi wavelength and I. 

Transport calculations have been able to explain the observed difference in 
the Weiss oscillatons for weakly ID and 2D modulated superlattices in a strong 
magnetic field. || The periodic potential lifts the degeneracy of the Landau 
levels (LL's) leading to Landau bands (LB's) of oscillating width periodic in 
J5 -1 . The bandwidth is determined [§|-|| by the commensurability of R c and 
L. 

The single-particle energy spectrum for this problem has been investigated 
by several researchers culminating in Hofstadter's butterfly, |l3|] a graph 

showing the complicated self-similar splitting of a LB into minibands as a 
function of the magnetic flux through a unit cell of the lattice. The effects of 
the splitting of LB's into minibands has, to the best of our knowledge, not 
been observed directly in the currently common superlattices (L ~ 300 nm 
n s = 3x 10 11 cm -2 ) since usually several LB's are occupied, which may overlap 
near Ep. Transport calculations based on the Hofstadter energy spectrum and 
taking into account short-range impurity scattering have, however, succeeded 
in recovering all the observed features of the magnetotransport coefficients for 
large-period superlattices. |J In the limit of a very high magnetic field the 
intricate gap structure of the Hofstadter butterfly is robust but long-range 
disorder tends to smear out the fine structure in the density of states. [15] 

The Hofstadter butterfly is obtained in two complimentary but mathe- 
matically equivalent limits: [[llj] first, in a strong lattice potential and a weak 
magnetic field in the tight-binding method, |3,11,13| and, second, a weak peri- 
odic perturbation of Landau quantized 2DEG. JD .14,16] Both approaches de- 
liver the so-called Harper's equation PJTl|] determining the energy spectrum. 
The energy spectrum has been investigated in the intermediate region in the 
absence of collision broadening. The coupling of the LL's by the external peri- 
odic potential strongly reduces the original high symmetry of the Hofstadter 
spectrum but retains a very complicated subband structure. [17,18] It has 
been possible to transform the equation of motion, the Schrodinger equation, 
into a vectorial form reminiscent of Harper's equation, but exact under the 
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most general conditions and exhibiting chaos in the classical limit. [18] 

All the above-mentioned theoretical investigations of the energy spectrum 
of a 2DEG in a periodic potential and a homogeneous perpendicular magnetic 
field have neglected the effects of the electron-electron interaction. In the 
present publication we initiate, within the Hartree approximation, the study 
of these effects on the electronic energy spectrum and the thermodynamic 
density of states. 



II. MODEL 

To describe the electrons in the conduction band of a lateral superlattice 
at the AlGaAs-GaAs interface in a constant perpendicular magnetic field B = 
Bz we employ a model of strictly two-dimensional electron gas (2DEG), with 
the three-dimensional charge density given by — en s (f)5(z), and f = (x,y). 
The orthogonal superlattice is spanned by the lattice vectors 

R = mli + nl2, n,m£Z. (1) 

l\ = l\x, and I2 = hy are the primitive translations of the Bravais lattice B. 
The corresponding reciprocal lattice 1Z is spanned by 

G = Gi5i + G 2 52, Gi,G 2 €Z, (2) 

where 

5i=2vr^, <? 2 = 2vrf- (3) 

The external periodic potential the electrons are moving in is taken to be of 
the simple form 

V(f) = V x cos(gix) + V y cos(g 2 y). (4) 

The electron-electron interaction is included in the Hartree approximation 
(HA) leading to the effective single-electron Hamiltonian 

H = H + Vh{t) + V(r), (5) 
where Vh(t) is the effective potential in a medium with a dielectric constant 

K, 

- r 2 , /Mr") - n h 



V H (f) = - / Sr' '"° , (6) 

K jb? \r — r 'I 

felt be each electron and caused by the total charge density of the 2DEG, 
— en s (r), and the neutralizing background charge density +enf, = +e(n s (r)). 

The periodic external potential V(r) and the constant external magnetic 
field imply that all physical quantities of the noninteracting system are peri- 
odic with respect to translations of R G B. Since the Coulomb interaction is 
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given by a central potential between pairs of electrons, the periodicity of the 
system is not broken by the electron-electron interaction. The periodicity of 
the Hartree potential @ follows from that of n s (r). 
The translation operator 

T(R) = exp(~R-ij, (7) 

where p is the momentum canonical to r, does not commute with Hq, 
[T(R),H ] / 0. Due to the presence of the vector potential A in H , 

Ho = ^-(p + -l) 2 , B = VxA (8) 
2m \ c J 

an additional gauge transformation is necessary. Defining the magnetotrans- 
lation 

S(R)=exp(^- X ^T(R), (9) 

where x is obtained from T _1 AT = A + Vx, we have for the mechanical 
momentum tt = p + (e/c)A and the position r 

S- l (R)rS(R) = r + R 

S-\il)nS(fi)=ir. (10) 

The Hamiltonian H is thus invariant with respect to the magnetotranslations, 
[S(R),H]=0. 

Hq defines the natural length and frequency scales for the system by the 
magnetic length / and the cyclotron frequency uj c 

eB 

l = \^ u c = —. (11) 

mc 

Generally the magnetotranslations (Q) do not commute, [S (Ri) , S (R2)] 7^ 0, 
expressing the fact that I is not commensurate with the lattice. But if an 
integer number of flux quanta $0 = (hc/e) flows through the lattice unit cell 
l\ x I2, then [S (R\) , S (R2)] = and the magnetic field assumes the values 

B = g l -J^ Y -, 9l eZ. (12) 
\h x h\ 

For these values of B the magnetotranslations (j9|) of the lattice vectors com- 
mute with each other, with H , and also with Hq. Thus, the known eigenvec- 
tors of Hq can be chosen to be also eigenvectors of the unitary magnetotrans- 
lations with simple eigenvalues 

S(t 2 )^ = e ie ^, Ote [-*,*]. (13) 
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The above eigenequations can be reinterpreted as periodic boundary condi- 
tions for the wave functions within one unit cell of B, and the eigenfunctions 
form a complete basis for each pair of phase angles (6±, 82)- 

Here we use the symmetric basis functions constructed by Ferrari [19] 
and used by Silberbauer. 20|, 2 1 ] A finer lattice with only one flux quantum 
through a unit cell is defined by the primitive vectors 

c = — , d = — , pq = gi,p,q£N. 
P q 

With the symmetric gauge A = (B/2)(— y,x) the functions 



-, 00 

(r) = — = E [S®* 

V m n=— 00 



S(d)e- l » <p ni (r), 



where 



<f>n, (r) 



6»i + 27rni 
P 

62 + 2im2 

q 
1 



y/2-irni\l 2 \ V21 



meh = {o,...p- 1} 

n 2 € h = {0, . . . q - 1} 



exp 



4|2 - n/ = 0,1,2..., 



form a complete orthogonal basis in the Hilbert space TCo^ if (/i, v) 7^ (vr, ir) 
for all (m, 71-2) £ I\ x The norm of the can be shown to be 



E (-i)™^ 



mn i(/im+wi) 



exp 



m,n=— 00 



1 

IP 



|nc + md|" 



and they have on the finer lattice the Bloch-type property 

In order to utilize the Ferrari-basis (|l5| ) {</> a }, each element of which satisfies 

H Q (j) a = E a (f> a , a = (/i, v, 111), 

to find the energy spectrum {e a } and the eigen functions {ip a } of H, 

H4> a = (H + V + V H )ip a = e a ip a , 

we need the matrix elements of the external periodic potential V(r) and the 
Hartree potential @ in that basis. Both V(r) and Vff(r) are periodic with 
respect to R £ B and can be expressed in terms of a Fourier series 

V{f) = E V(G)e^, 
Gen 



with the Fourier coefficients given by 



V(G) = - - f d 2 r V(r)e- iT "- d . 
A J a 

Here A stands both for the unit cell in B and it's area \l± x Z2 1 . The task 
of deriving the matrix elements of the external potentials is thus reduced to 
finding the matrix element of the Fourier phase factor. This has been accom- 
plished by Silberbauer showing that the matrix elements can be expressed 
as rapidly converging series of simple analytic functions. [20,^1]] The result is 



reproduced here, correcting a simple typographical error in one of the sources 
fl2C 



Kn^|e iG >in 2 n z ) = ^(G)TS(G)-^ 

\\<0 2 ||||0nH 
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where 



4/ 2 



T$$(G)= £ (-l)^ e ^'A+^)exp|-^G(A C + ^r l Ar + "-' / l" 

AQ=-oo I 

and 

G n 'i ni {G) -- 



\G\ z r 



ni 



Here a complex notation, G = G\\gi\ + iG^I^I etc., has been used for the 
vectors G, Eq. and similarly for c, and d, Eq. are the associate 

Laguerre polynomials. Furthermore, the matrix elements T n *ni(G) vanish 
unless G\ = rij — n\ — Mp and G2 = n' 2 — ri2 — Nq for some M, N 6 Z. 
The Fourier coefficients for the Hartree potential are 
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, 2ir^-An s (G) if G^O 
V H (G) = { 1 ; 

if G = 



where 



1 



2„/ fz*t\ „ } „-ir'-G 



An s {G) = - I d 2 r' {n s (r')-n b } 



e 



A 



The charge neutrality of the system removes the singularity of the Coulomb 
interaction for G = 0. 

To complete the Hartree set of equations (|j) , (||) , (g) , and (||) the electron 
density is written in terms of {e a } and {ip a } as 



n *( f ) = ^fE£ f(e a - fi)\Mr)\' 



where / is the equilibrium Fermi distribution and \x is the chemical potential. 
This set of Hartree equations has to be solved iteratively together with the 
condition that the average electron density n s = N s /A is constant and deter- 
mines the chemical potential fj,. N s is the number of electrons per unit cell 
with area A in B. 
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III. RESULTS 



The Coulomb interaction couples directly the subbands of a particular 
LL and the subbands of different LL's via the matrix elements (|23|). The 
matrix elements with states of different phase angles (9%, 62) do all vanish, 
so that the Hartree Hamiltonian can be diagonalized for each pair (9±, 62) 
separately. However, the states for a given set of phase angles depend on the 
states of all other phase angles via the Hartree potential (|6|) and the density 
n s (r) (|28|). The self-consistent model described in the above section can be 
used to evaluate the energy spectra {e a } and the wavefunctions {ip a } for a 
wide variety of modulation strengths V x and V y and values of the magnetic 
field B by selecting the appropriate size of the basis, H is diagonalized in. If 
only the basis states in the same LL are taken into account, the modulation 
strength is irrelevant since it can be factored out of the Hamiltonian matrix. In 
this situation, corresponding to the usual discussion of Hofstadter's butterfly 
on the basis of Harper's equation, a highly symmetric energy spectrum is 
expected. Of course, this restricted model is only appropriate for describing a 
system with very weak modulation (as compared with the inter-Landau-level 
energy spacing huj c ). In the following, parameters appropriate for GaAs have 
been used, i.e. m* = 0.067m-o, k = 12.4, and the spin degree of freedom 
has been completely neglected since the exchange force is not accounted for 
in the interaction between the electrons. The superlattice is selected to have 
square symmetry, l\ = h = L, and V x = V y = Vo, to comply with the 
usual experimental conditions, even though the Ferrari-basis is applicable for 
a much more complicated unit cell. [|lS|-|2l|l 

The iteration of the Hartree equations (||), (|5|), (||), and (g) has been 
accomplished by specifying an initial electron density that is then used to find 
the Hartree potential (0). Subsequently the Hamiltonian (]5|) is diagonalized 
yielding {e a } and that are in turn used to derive a new density (|28|) . In 
order to guarantee convergence the speed of the iteration process is reduced 
by using only 10% of the new densty and 90% of the density from the former 
iteraction cycle. The convergence is considered completed when the integrated 
rms deviations of the density from two subsequent cycles falls below 0.01% 
of N s . Clearly the convergence rate depends on modulation stength of the 
external potential (||), the filling factor v and the temperature T. Commonly 
for the following parameters we need between 15 and 30 cycles to attain 
the specified accuracy. For the ground state properties the evaluation of the 
wave functions can be avoided by combining equations ( |27|) and (p8| ) and thus 
write the Fourier transformed density entirely as a function of the analytically 
known matrix elements (|23|). 

The energy spectrum for the case of four flux quanta $0 piercing the unit 
cell of a Hartree interacting 2DEG periodically modulated with L = 50 nm 
is shown in Fig. The dimensionless magnetic flux is thus, <&/<&o = pq = 4, 
which translates into a magnetic field of B ~ 6.617 T and huj c ~ 11.43 meV. 
Due to the strong modulation of the periodic potential (||), V$ = 4.0 meV, 
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and the low temperature T = 0.2 K, four separate subbands of the lowest LB 
are resolved, which are centered about E = fiui c /2 ~ 5.7 meV. The electronic 
density is very low, N s = 0.25 being the average number of electrons per unit 
cell. The chemical potential /x = 2.38 meV is lying in the lowest subband 
which has a very slight dispersion as function of the phase angles 9\ and 
#2- This is a manifestation of the familiar pinning effect observed earlier in 
many different systems. [22],23| The relatively weak screening (judged from 
the spread of the minibands < 2Vo) is strongest in the subband where \x is 
situated. The next LB is centered around E = 3Huj c /2 ~ 17.1 meV and 
has a much weaker dispersion than the first one. The vanishing overlap of 
the LB's leads to an approximate electron-hole symmetry between the cases 
N s = 0.25 and N s = 3.75, since the (average) filling factor v = 1 is reached 
for N s = 4. Note however that, e.g. for N s = 0.25, the energy spectrum 
is not symmetric around the mean energy, in contrast to the case of the 
noninteracting electrons. This approximate symmetry can also be seen in 
Fig. |2] in which the self-consistent potential V sc (f) = Vuif) + V{f) for pq = 6 
has been depicted. Here the intersection of the shifted chemical potential fi — 
huj c /2 with V sc projected onto the basis of the graph is also shown. For N s = 
0.25 (Fig. |a) the 2DEG forms isolated lakes near the minima of V sc . The 
complimentary case, N s = 5.75, is depicted in Fig. |2|c, where the empty states 
of the first LB, "the holes," are isolated near the maxima of the potential. The 
transition between these two regions is seen in Fig. ^o, where traditionally, i.e. 
in a semi classical picture, neither the electrons nor the holes are thought to 
form closed orbits. The screening is so strong for this last case (corresponding 
to v = 1/2) that V sc is almost flat. As stated earlier, the symmetry is only 
approximate, since the states of the second LB are mixed into the states of 
the first one by the interaction. The potentials are thus not identical for 
the complimentary cases of N s = 0.25 and 5.75, and the transition between 
the two regions does not happen exactly at N s = 3.00. The self-consistent 
densities of the electrons for the cases N s = 0.25 and N s = 5.75 is shown in 
Fig. |3[ underlining clearly the near isolation of the electrons in the minima of 
V sc when N s = 0.25. In Fig. || we compare the self-consistent density n s {r) 
with the density of the non interacting 2DEG for N s = 3.00. For pq = 6 
the screening is at its maximum strength but due to the strong modulation 
Vo = 4.0 meV the interaction only causes a moderate flattening of the density. 

The strong screening encountered in Fig. |2| for the pq = 6 case is even more 
drastic for a superlattice potential V(f) with a longer period and weaker 
modulation as can be observed from Fig. |5| when L = 200 nm and Vq = 
0.5 meV. Here the width and the location of each subband of the lowest 
LB is shown as a function of N s , the location of the chemical potential is 
also indicated. Even though the temperature is quite low, T = 0.1 K, the 
Coulomb interaction between the electrons is very effective in screening the 
superlattice potential and collapsing the subbands into an almost degenerate 
LL around N s ~ 3 corresponding to v ~ 1/2. For an integer filling factor 
of the order of unity (N s ~ 6) the screening is as expected very weak. [22] 
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The upper panel of Fig. [B| shows the path of the chemical potential n through 
the subbands as a function of the average filling factor u or N s when the 
electrons do not interact. The fine structure of \x leads to a peak structure of 
the (dimensionlessly written) thermodynamic density of states (TDOS) 

-w 1 , n s = (n s (r)) 

whenever fi crosses a subband. The strong screening in the case of the inter- 
acting 2DEG (the lower panel of Fig. ||) results in only one peak in Dt, in 
which no fine structure is resolved. The corresponding information is shown 
for four different cases in Fig. ||, pq = 2, 3, 4, 6, identifiable from the number 
of subbands appearing in each figure. Here the period is short, L = 50 nm, 
the modulation strong, Vo = 4.0 meV, and T = 1.0 K. However only one LL 
is used in the calculation, so that the results apply to the case of a weak su- 
per lattice potential. Here the electron- hole symmetry is exact, as is reflected 
in each of the four cases by the point symmetry of the subband structure 
around the chemical potential /i at half filling. The screening is strongly de- 
pendent on the total number of states in a lattice unit cell and the average 
filling factor. It is therefore, strongest for pq = 6 when maximum six electrons 
can occupy the lowest Landau level in the unit cell and when the LB is half 
filled, N s = 3. The thermodynamic density of states Dt, Eq. (p9|), for the 
corresponding cases is shown in Fig. |7]. Clearly the strong screening in the 
pq = 6 case obscures the fine structure in Dt leading to only one pronounced 
peak. The subbands in the other cases turn up as separate peaks in Dt- The 
calculation has been repeated for the same parameters but one more Landau 
level resulting in qualitatively the same picture. The electron-hole symmetry 
is then only approximate with the largest deviation for pq = 2, 3, but the 
resolution of the peaks in Dt is unchanged. 

The Hofstadter butterfly should be recovered when the calculation is per- 
formed for N s = and only one LL, and if the width and location of the 
subbands are plotted as functions of the inverse of the number of magnetic 
flux quanta through a unit cell of the lattice, i.e. of ^o/^ = ^/(pq)- Ln 
addition, all energies have to be scaled according to: E — > (E — fru c /2)/Vo. 
For interacting electrons the number of electrons per unit cell N s , or the 
average filling factor of the LB turns out to be a new parameter control- 
ling the bandstructure. We therefore, present in Fig. || the subband struc- 
ture for the four inverse fluxes l/(pq) = 1/2,1/3,1/4,1/6 in the four cases 
N s = 0.00,0.25,1.00,1.50 together with the location of the chemical poten- 
tial fi (except for N s = 0.00, where it is irrelevant). For low density, say 
N s = 0.25, the subband structure is almost symmetric around the energy 
zero, like in the noninteracting case, but the subbands become quite asym- 
metric for a higher density of electrons. For the short-period super lattice 
studied here, and for a low density of electrons, the essential gap structure 
does survive in the presence of interaction. The presentation in each of the 
subfigures of Fig. [8| corresponds to the experimental procedure of keeping 
the density of electrons fixed but changing the magnetic field. Another way 
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to investigate the screening is to keep the average filling factor v constant 
but vary the magnetic field and, thus, also the average density of the 2DEG. 
Fig. |9| compares the subband structures for l/(pq) = 1/2,1/3,1/4,1/6 and 
v = 1/2 with the complete Hofstadter butterfly. Here the energies have been 
scaled with the factor (E — hu c /2) exjp{(-irl / L) 2 } /Vo so that the results can 
be directly compared with an earlier calculation of the Hofstadter spectrum. 
The energy spectrum in the interacting case shows an overall reduction 
in dispersion or bandwidth due to the strong screening that is most effective 
for large flux and large number of available states. The bandwidths of the 
subbands for the interacting case has been evaluated here on a discrete lattice 
in the (6\, #2)-plane without attempting an interpolation between the lattice 
points; thus, the actual bandwidths can be larger by a small percentage of the 
widths shown. This effect explains why the bands for the even denominators 
in Fig. B for N s = are not touching. 



IV. SUMMARY 

In an interacting 2DEG subject to a super lattice potential and a homoge- 
neous perpendicular magnetic field not only the magnetic flux through a unit 
cell but also the density of electrons determines the complicated splitting of 
the Landau bands into subbands. We have shown that in the Hartree approx- 
imation the essential gap structure of the energy spectrum remains, although 
the screening leads to a quenching of the Hofstadter butterfly at small values 
of the inverse flux (see Fig. 8). The symmetry of the butterfly is also low- 
ered by the coupling to higher LL's due to increased strength of the periodic 
potential, as has already been discussed by other authors. []T7| , p^ 1 For periods 
around L = 200 nm (currently attainable in experiments on super lattices) the 
2DEG can effectively screen the periodic potential even for a very low den- 
sity n s so that only one or two Landau bands are partly occupied. Only at 
shorter lattice constants (L < 100 nm) and thus much higher magnetic fields 
we can, on the basis of our Hartree calculation, expect the subband structure 
to be resolvable in experiments, when n s is maintained low enough. These 
predictions are made on the basis of the calculated energy spectrum and the 
structure observed in the thermodynamic density of states Dt, see Eq. (|29|). 
The transport coefficients do not depend on Dt in any simple way, and may 
be more sensitive to the subband structure of the energy spectrum. M In 
addition, the collision broadening due to impurities and inhomogeneities has 



to be considered, especially, when the electronic density is very low. |3,24,15| 
It is to be expected that the inclusion of exchange and correlation effects 
may lead to more dramatic changes in the energy spectrum than the direct 
interaction considered in the Hartree approximation. One might even expect 



rearrangement of the bands or the occurence of spin-density waves. [ 25, 23, 26 1 
A reliable procedure for the treatment of exchange and correlation effects in 
the 2DEG in a strong magnetic field and a lateral superlattice is, however, at 
present not available. Therefore we have restricted our investigation of the 
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electron-electron interaction effects by considering the Hartree approximation, 
which is expected to yield a qualitatively correct description of the dominant 
screening effects in inhomogeneous electron systems. 

In order to keep the computational efforts in reasonable limits, we have 
restricted our calculations to a few characteristic values of the magnetic field. 
Our results indicate in which manner screening effects will change the over- 
all appearence of the Hofstadter energy spectrum. Unfortunately, there is 
no simple way to extrapolate from these special values to arbitrary rational 
numbers of flux quanta per unit cell. For a more detailed comparison with 
the energy spectrum of the noninteracting case, further studies using a much 
denser set of flux values are desirable. However, the answer of the interest- 
ing question, whether the self-similarity of the Hofstadter butterfly survives 
in the presence of the mutual Coulomb interactions between the electrons, 
at least within the Hartree approximation, demands an amount of computa- 
tional work which seems at present unfeasible. On the other hand, from the 
experimental point of view these details of the energy spectrum seem not to 
be accessible anyway. The challenge for the near future is to resolve the most 
prominent gaps of the energy spectrum experimentally. 
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FIGURES 



FIG. 1. The dispersion of the four minibands of the lowest Landau level in the energy spectrum 
of the Hartree interacting 2DEG as a function of the phases 9\ and 62 when four flux quanta of the 
magnetic field pierce a rectangular unit lattice cell (p = q = 2). L x = L y = L = 50 nm, T = 0.2 
K, V x = V y = Vq = 4.0 meV, N s = 0.25, and fx = 2.38 meV. Two Landau levels are included in the 
calculation. 

FIG. 2. The self-consistent potential V sc (r) = Vh(t) + V(r) in a unit lattice cell for (a) 
N s = 0.25, (b) N s = 3.00, and (c) N s = 5.75. The intersection of the chemical potential \i — hu c /2 
and V sc is projected on the base and [i — hu> c /2 is shown in the upper right corner of each subfigure 
in meV. p = 2, q = 3, L = 50 nm, T = 1.0 K, Vq = 4.0 meV. Two Landau levels are included in the 
calculation. 

FIG. 3. The electron density n s (r) (nm' 2 ) in a unit lattice cell for (a) N s = 0.25, and (b) 
N s = 5.75. Contours of n s (f) are shown at the base, p = 2, q = 3, L = 50 nm, T = 1.0 K, 
Vq = 4.0 meV. Two Landau levels are included in the calculation. 

FIG. 4. The electron density n s (r) (nm' 2 ) in a unit lattice cell for N s = 3.00. (a) Non 
interacting, and (b) Hartree interacting 2DEG. Contours of n s (r) are shown at the base. Same 
scales are used in both cases, p = 2, q = 3, L = 50 nm, T = 1.0 K, Vq = 4.0 meV. Two Landau 
levels are included in the calculation. 

FIG. 5. The bandwidth of the minibands of the lowest Landau band as a function of the 
number of electrons per unit cell N s for a noninteracting 2DEG (upper) and a Hartree interacting 
2DEG (lower). The chemical potential ji is noted by a continuous curve. L = 200 nm, T = 0.1 K, 
Vq = 0.5 meV, p = 2, and q = 3. Two Landau levels are included in the calculation. 

FIG. 6. The bandwidth of the minibands of the lowest Landau band as a function of the number 
of electrons per unit cell N s for a Hartree interacting 2DEG and (p = 1, q = 2) (upper left), (p = 1, 
q = 3) (upper right), (p = 2, q = 2) (lower left), and (p = 2, q = 3) (lower right). The chemical 
potential u is noted by a continuous curve. L = 50 nm, T = 1.0 K, Vq = 4.0 meV. One Landau 
level is included in the calculation. 

FIG. 7. The thermodynamic density of states l 2 hu> c (dn s /dfi)TB as a function of the dimension- 
less inverse magnetic flux through a unit cell l/pq and N s . L = 50 nm, T = 1.0 K, Vq = 4.0 meV. 
One Landau level is included in the Hartree calculation. 

FIG. 8. The scaled bandwidth of the subbands of the lowest Landau band (E — Uui c /2) /Vq as 
function of the dimensionless inverse magnetic flux l/pq for N s = 0.00 (upper left), N s = 0.25 
(upper right), N s = 1.00 (lower left), and N s = 1.50 (lower right). The chemical potential u is 
noted by crosses. L = 50nm, T = 1.0K, Vq = 4.0meV. One Landau level is included in the Hartree 
calculation. 
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FIG. 9. The scaled bandwidth (E — fwj c /2) exp{(7r//L) 2 }/Vo of the subbands of the lowest 
Landau band as function of the dimensionless inverse magnetic flux 1 /pq for the Hartree interacting 
2DEG with v = 1/2 (left), and for the noninteracting 2DEG (right), the Hofstadter butterfly. One 
Landau level is included in the Hartree calculation. In the left subfigure the 6 LB's for the case 
pq = 6 can not all be resolved due to vanishing band gaps. The parameters used for the left 
subfigure are: L = 50 nm, T = 1.0 K, Vq = 4.0 meV. 
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